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ABSTRACT 

It has been shown that in an idealized galaxy cluster survey, containing 10, 000 clusters, statisti- 
cal errors on dark energy and other cosmological parameters will be at the percent level. Furthermore, 
through "self-calibration" , parameters describing the mass-observable relation and cosmology can be 
simultaneously determined, though at a loss in accuracy by about an order of magnitude. Here we 
examine an alternative approach to self-calibration, in which a parametrized ab-initio physical model 
is used to compute theoretical mass-observable relations from the cluster structure. As an example, 
we use a modified-entropy ("preheating") model of the intracluster medium, with the history and 
magnitude of entropy injection as unknown input parameters. Using a Fisher matrix approach, we 
evaluate the expected simultaneous statistical errors on cosmological and cluster model parameters. 
We find that compared to a phcnomcnological paramctrization of the mass-obscrvablc relation, our 
physical model yields significantly tighter constraints in both surveys, and offers substantially im- 
proved synergy when the two surveys are combined. In a mock X-ray survey, we find statistical errors 
on the dark energy equation of state are a factor of two tighter than the phenomenological model , with 
Awq ~ 0.08 and its evolution, Aw a = —Adw/da ~ 0.23, with corresponding errors of Awq ~ 0.06 
and Aw a ~ 0.17 from a mock Sunyaev-Zel'dovich (SZ) survey, both with N c i ~ 2.2 x 10 4 clusters, 
while simultaneously constraining cluster model parameters to <J 10%. When the two surveys are 
combined, the constraints tighten to Awq ~ 0.03 and Aw a ~ 0.1; a 40% improvement over adding the 
individual experiment errors in quadrature, and a factor of 2 improvement over the phenomenologi- 
cal model. This suggests that parametrized physical models of cluster structure will be useful when 
extracting cosmological constraints from SZ and X-ray cluster surveys. 

Subject headings: cosmology:observations-cosmology:theory-cosmological parameters-dark matter- 
galaxies:clusters:general-large-scale structure of the universe 



1. INTRODUCTION 

It has been known for about two decades that clus- 
ters of galaxies provide a uniquely powerful probe of 
fundamental cosmological parameters. As the most 
massive virialized objects in the universe, they form at 
the highest peaks in the primordial density field, and as 
a result, their abundance and spatial distribution are 



very s ensitive to the underlying cosmology (e.g.JBahcal] 
1988; Pe ebles! H99l ICarlberg et all H997t [ 



2002; 



Peebles! 119931: ICarlberg et~ai1 Il997t iRosati et al 
Henrvi 20041 ). Strong constraints have already 



been derived on the present day matter density O m , and 
the normalization of the dark matter p ower spectrum as 
(iPeebles. Daly, fc Juszkiewiczl 1 19891 iHenrv fc Arnaudl 
19911: iBahcall fc Cenl fl992: lyTanafcLiddlel il99ft 



Bahcall fc Bode! l2003t IBahcall et afl l2003[ h using only 
dozens of clusters. In large future surveys, with tens 
of thousands of clusters, high-precision constraints are 
expected to be available on d ark energy parameters 
([Haiman. Mohr. fc Holder! [2001). including constraints 
on the evolution of the equation of state w a = dw/da 
(IWeller. Battve. fc Kneissl l2002t IWeller fc Batty3 l2003t 
iWang et al.ll2004l ). 

Plans are are currently being drafted for the first very 
large scale clusters surveys. These projects, such as the 
4,000 degree 2 South Pole Telescope that will find clusters 
through the SZ effect (hereafter SPT; Ruhl et al. 2003), 



a proposed larger X-ray survey (| Haiman et al. I [20051 ). or 
catalogs of clusters i dentified in w ide weak-lensing sur- 
veys, such as LSST ([Tvsonl [20021 ). promise to produce 
catalogs of many thousands of clusters, containing a rich 
resource of cosmological information. With such an in- 
crease in the data looming on the horizon, it is becoming 
increasingly important to understand the inherent sys- 
tematics, and establish realistic expectations for poten- 
tial cosmological constraints. 

Several recent works, motivated by the above consid- 
erations, have focused on various aspects of extracting 
cosmological parameters from high-yield future surveys, 
such as the statistical constraints avail able on curva- 
ture fife ([Holder. Haiman. fc Mohri |2001|). assessin g the 
impact of sample v ariance dHu fc Kravtsovl l2003f) and 
other uncertainties ([Levine. Schultz. fc White! l2002h on 
parameter estimates, and controlling such uncertainties 
by utilizing information fro m the sh ape of the clus- 
ter mass function dN/dM ((Hu| l2003f ). Recent stud- 
ies have also elucidated the additional cosmological 
information available from the spatial distribution of 
galaxy clusters through a me asurement of their t hree- 
dimensional power spectrum ([Hu fc Haiman! I2003T ) . uti- 
lizing _both^_the_jntr^ func- 
tion (iRefregier. Valtchanov. fc Pierre! I2002D and baryon 
features dBlake fc Glazebrookj 120031 : ISeo fc Eisensteinl 
l2l)03 Hinder 200|). 
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Most importantly, recent work has focused on the sys- 
tematic errors on cosmological parameters arising from 
the inherent uncertainties on the mass-observable rela- 
tion. It has been shown that large forthcoming cluster 
surveys can be "self-calibrating" , in the sense that when 
the abundance and clustering properties of clusters are 
considered in tandem, tight constraints can be derived 
on cosmology even if the mass-observable relation has 
to be determined simultaneously from the same data 
(iMaiumdar fc Mohri[200l iWang et alMooi [Lh^a~fcHul 
|2005| ). In these studies, the mass-observable relation was 
assumed to be either unknown apriori, or to have a sim- 
ple parametric form (such as a power-law). 

While these assumptions have been demonstrated to 
yield interesting constraints, one may argue that this ap- 
proach is overly conservative. In particular, it is rea- 
sonable to expect that the structure of clusters, such 
as their mean density and temperature profiles, can be 
at least approximately computed from ab-initio models. 
The simplest, spherically symmetric, self-similar mod- 
els contradict observations, such as the relation between 
X-ray flux (L x ) and tempera ture (T), or the resolved 
profiles of low-mass clusters (jVoit et al.l l2002D. How- 
ever, motivated by the observations of iPonman et al.l 
(1999), modified-entropy models of the intra-cluster 
medium (ICM), whi l e still not correct in detail 
(iPratt fc Arnau d 2005; Pratt, Arnaud. fe Pointecouteaul 
20061; lYounger fc Brvanl 120061 ) . have been successful at 
predicting the mass-observable scalings in te mperature, 
luminosity, and SZ effect, from first principles (jVoit et al.l 
l2002t iMcCarthv et al.ll2003al lbl). It is reasonable to ex- 
pect that our understanding of the physics determining 
the structure of clusters will be improved in the future 
by advances in both observations and numerical simula- 
tions. 

In this paper, we study the utility of a partial un- 
derstanding of the ICM in improving cosmological con- 
straints. We will show that using a-priori cluster struc- 
ture models, it is possible to significantly improve cosmo- 
logical constraints, relative to direct, phenomenological 
self-calibration of the mass-observable relation - even if 
the physical model has free parameters. Furthermore, 
the combination of X-ray and SZ datasets is pa rticu- 
larly promising (|Verde. Haiman. fe Spergel |2002|) , and 
can more effectively break degeneracies, since the un- 
derlying models directly relate the SZ decrement to the 
X-ray flux. The end result is significantly improved er- 
rors on the cosmology, using the same observational data 
set, in addition to a well constrained model of the ICM. 

The rest of this paper is organized as follows. In 
we describe our models of cluster structure and evolu- 
tion in the presence of a preheating of the ICM gas. 
In Jj3l we discuss our fiducial model, and demonstrate 
that it satisfies a number of existing observations. In 
SI we contrast our physical models with a "traditional" 
conservative approach to self-calibration, in which the 
mass-observable relation is parametrized directly as a 
power-law. In we discuss and motivate the param- 
eters of our mock cluster surveys, and, for comparison, 
discuss future measurements of cosmic microwave back- 
ground (CMB) anisotropies. In we briefly summarize 
the Fisher matrix approach, as applied to cluster sur- 
veys and CMB anisotropies. In $3 we present and crit- 
ically discuss our results, by comparing the constraints 



obtained under self-calibration and under the assump- 
tion of physical cluster models. Finally, in HWl we sum- 
marize our main conclusions and the implications of this 
work. 

2. CLUSTER STRUCTURE AND EVOLUTION MODELS 

In this section, for completeness and specificity, we de- 
scribe our models for the structure - density and tem- 
perature profiles - of galaxy clusters, in the presence or 
absence of excess entropy. Our treatment assumes that 
the gas is in hydrostatic equilibrium in a dark matter 
(DM) potential, and it closely follows previous studies 
on cluster profiles and preheating. 

2.1. Motivation 

Modified entropy models are physicall y motivated by 
the pr esence of a universal entropy floor. IPonman et al.l 
(1999) find that the entropy distribution of ICM gas, 
as inferred from ROSAT observations of 25 galaxy clus- 
ters, is inconsistent with simple gravitational collapse, 
and suggest instead the presence of an entropy floor of 
Kq ~ lOO/i -1 / 3 keV cm 2 . This excess entropy is itself 
potentially a relic of winds and outflows from galaxies 
that formed in the protocluster, and heated the ICM at 
moderate redshift z ~ 2, prior to cluster collapse. They 
then argue that this results in the observed departure of 
real galaxy clusters from the expected self-similar scal- 
ing. However, it is plausible that only a fraction of this 
preheating entropy is injected at high redshift. Some of 
the heat input may come from a combination of active 
galactic nuclei (AGN) and supernovae (SN), over an ex- 
tended period, rather than winds in a narrow redshift 
range. For this reason, it is useful to treat Kq as a func- 
tion of redshift, which can be parametrized, for example, 
by a power-law relation K (z) = K (z = 0)(1 + z) aK . 

We assume that clusters at a given redshift z are spher- 
ically symmetric and in hydrostatic equilibrium, with an 
ideal gas ICM consisting of a fully ionized H-He plasma 
with the cosmic helium mass fraction Yhc = 0.25, and 
a mean molecular weight /.t = 0.59. The unmodified en- 
tropy distribution is constructed by assuming that the 
gas traces the dark matter (DM), and that the total 
gas density is given by a fraction f g of the cosmologi- 
cal baryon density within the virialized cluster that re- 
mains in t he ICM. The DM and gas profile is assumed to 
follow the lNavarro. Frenk fc White] (fl997h (NFW) form, 
as motivated by N-Body si mulations, with a fixed con- 
centr ation parameter c = 5 (Ek e. Navarro, fc Steinmetzl 
2001). The assumption of hydrostatic equilibrium, to- 
gether with physically motivated boundary conditions, 
allows us to calculate the pressure and temperature pro- 
files, which, combined with the gas density, give the en- 
tropy distribution as a function of radius. This entropy 
distribution is then modified by adding a constant Kq at 
all radii, and the hydrostatic equilibrium equations are 
re-integrated under the assumption that mass shells do 
not cross (i.e. K{M g ) is conserved where M g is the total 
gas mass at that radius). The resulting modified density 
and temperature profiles have flat cores, and they can 
be used to predict cluster observables, such as the X-ray 
flux or SZ decrement. 

2.2. Unmodified Equilibrium Structure 
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The ICM is assumed to trace the DM, which is well- 
described by an NFW profile of the form 

0~cPc 



Pdm{t) = 



(1) 



(r/R s )(l+r/R s )i 

where R s is a scale radius defined in terms of the con- 
centration paramet ers c and the virial radius, as in 
iBullock et ahl (l200l . by 



cR s = R v 



( 3M 



\Aml\ v p c 



1/3 



and 



Sr.= 



3 ln(l 



(2) 



(3) 



l+c 



In both cases A„ is the virial overdensity relative to the 
critical density of the universe p c given as a fitting form 
for a flat universe by iKuhlen et alj (|2005| ), extended to 
include a time-evolving dark energy: 



187r 2 fi m (z)[l + aO b (z)} 



(4) 



where Q(z) = Vl^iz) — 1, £l m {z) is the matter density 
of the universe as a function of redshift z, a = 0.432 — 
2.001(Hz)| - 234 - 1), b = 0.929 - 0.222(Hz)| - 727 - 1), 
and w comes from the dark energy ecjuation of state w = 
wo+Waj^ dChevallier fc Polarski|[200it iLinder! 120031) . 
The gas is assumed to trace the DM distribution for a 
given cosmology with a constant gas fraction /„ ~ 0.8, 



Pg( r ) 



ScPc 



(r/R s )(l + r/R s )' 



-fgfb 



(5) 



where = fit,/f2 TO is the mean baryon fraction of the 
universe. 

The corresponding ICM pressure and temperature pro- 
files are computed by integrating the equations of hydro- 
static equilibrium and mass conservation for an ideal gas, 



dP 



dM g _ 
dr 

P = nk B T 



4irr p g 

_ p g k B T 
pm p 



(6) 



(7) 



(8) 



where ks is the Boltzmann constant, m p is the mass of 
a proton, p = 0.59 is the mean molecular weight, and 



7(0 = - 



G(M 



DM 



Mg) 



(9) 



is the gravitational acceleration including the self-gravity 
of the gas. Once P(r) and T{r) are determined, the 
entropy profile K(r) is computed according to 



K = Tn- 2/3 



(10) 



where n is the total gas 1 particle number density n — 
P/ksT. The outer boundary condition is taken to be 
a pressure resisting the infall of baryonic matter at the 
virial radius; P„ir = fbpDMVft/3, with vjf the free-fall 
velocity at the virial radius as in IVoit et alJ (|200l . 

This in contrast to typical observational definitions which use 
the electron density n e 



2.3. The Modified Distribution 

After the unmodi fied entropy distr ibution has been cal- 
culated, we follow IVoit et al . (2002), and increase it by 
a constant Kq as 



K(r) = K(r) + K 



(11) 



where K(r) is the unmodified distribution and Kq is a 
constant of order 100 keV cm 2 , assumed here to evolve 
with redshift as Kq = Kq(z = 0)(1 + z)° K . This mimics 
the preheating of the ICM, whereby some mechanism had 
injected a fixed amount of entropy into the gas prior to 
its co llapse at redshift z (|Ponman et al.|[l999l : IVoit et al.l 
[2002]) . 

We change variables and re-integrate the equations of 
hydrostatic equilibrium using the modified entropy dis- 
tribution K. The equation of hydrostatic equilibrium 
then reads 



dP 

dr 



= g(r)p g (P,K) 



(12) 



where, casting the pressure and temperature in terms of 
P and K 



p g = pm p 



p 



-1 3/5 



k B K(M g 



(13) 



and 



k B T = K{Mgf /b P 2/b 



(14) 

We use the same pressure matching boundary condi- 
tion as before. This procedure results in modified profiles 
for P, pg, and T that can be used to make predictions for 
different bulk properties of the cluster, and to establish 
the expected mass-observable relations. 

2.4. Radial Profiles 

Before considering the effects of entropy modification 
on the bulk properties of the ICM, it is instructive to 
visualize its impact on the radial profiles of clusters of 
different masses. In particular, for reference, we examine 
the effects of entropy modifications with Kq = 50 and 100 
keV cm 2 as compared to the unmodified model on two 
clusters with M V i r — 10 14 and 10 15 h^ 1 M & respectively, 
at redshift = 0, with our fiducial ACDM cosmology 
(see definition below). The radial temperature, pressure, 
and gas density (T, P, p g ) profiles are plotted in Figure[TJ 

The overall impact of the entropy injection is to depress 
and flatten the core pressure and densi ty. This results 
in a flattened inner gas d ensity profile (|Sun et alJl2003t 
iPointecouteau et al 2004) and an outer (r > 0.1i?„) tem- 
perature profile ( Vikhlinin et al.ll2004h . resembling the 
results of recent X-ray observations. It is important to 
note that the effects of preheating are more pronounced 
for a low mass cluster than for a high mass cluster, with 
the limit that at very high mass the entropy injection 
becomes unimportant. This is because the injection is 
a decreasing percentage of the unmodified entropy at a 
given radius as the cluster mass increases. 

2.5. Observables 

In order to judge the fitness of our model for use in 
later cosmological measurements, we first need to confirm 
that it can reproduce the observed scalings for the M — 
T and L x — T relations, given the right choices for the 
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input parameters f g and Kq. We recall that the X-ray 
luminosity scales approximately as L oc p 2 , and define 
the luminosity weighted temperature as 



p 2 Jr)T(r)dV 



p 2 Ar)dV 



(15) 



where the integrations (as in all the following cases) are 
performed over the entire volume of the cluster out to 
the virial radius r — R v . 

We then assume a fully ionized H-He plasma with a 
hydrogen mass fraction X = 0.75, n e — Pg{^~ + ir^T") 

and nu — Pg^~i an d calculate the X-ray luminosity as 



L, 



n e riH A„(T)dV dv : 



(16) 



Here A„(T) is the iRavmond fc Smithl ()1977f ) cooling 
function for a gas of metalicity Z = 0.3Z Q for either 
the bolometric or the K-corrected soft X-ray emissivity 
in the 0.5 - 2.0 keV band. 

Inverse Compton scattering by electrons in the hot 
ICM manifests itself as a fractional change in the tem- 
perature of the Cosmic Microwave Background (CMB), 
which at fre quencies less than v ~ 218 GHz appears as 
a decrement (jSunvaev fc Zeldovich|[T97l I1980D . This ef- 
fect, often referred to as the thermal SZ effect, provides a 
useful probe of the electron pressure P e integrated along 
the line of sight (assuming the electrons and gas are in 
thermal equilibrium), which complements the X-ray lu- 
minosity measurements. We calculate both the central 
and integrated SZ flux according to the prescription in 
IMcCarthv et alj (|2003aUbl ). 

The amplitude of the SZ effect is encoded in the Comp- 
ton parameter, 



y(0) 



m e c 2 



P e dl 



(17) 



where <jt is the Thompson cross section and P e — n e kBT 
is the electron pressure integrated along the line of sight 
dl. The central decrement yo = y(9 = 0), and the SZ 
flux is 

Sv=jv{x)y in t (18) 

where x = hv I UbTcmr with Tqmb = 2.728 
(|Fixsen et alj I1996D is the dimensionless frequency, 
j v {x) — 2(kBTcMB) 3 {hc)~ 2 f v (x) describes the spectrum 
of the SZ effect, and 



Ifent(< 9) = 2tt / y(0')0'd0' 
Jo 



(19) 



We choose our observable as the integrated SZ flux, 

2(k B T C MB) 3 , , on x 
Su = Vint jj-^ U, (20) 

evaluated at the frequency v — 150 GHz. 

3. THE FIDUCIAL MODEL AND COSMOLOGY 

In this section, we discuss the parameter choices in 
our fiducial model. Our aim is to demonstrate that our 
fiducial model is in rough agreement with a number of 
existing observations. 

The observations we are mainly interested in match- 
ing are the scaling relations between cluster mass and 



observables, such as M^qq — L^oi, -M500 — S v Jf v , and 
Lboi — Su/fis, where Ma represents the mass within a 
spherical overdensity of Ap c . These relations drive the 
constraints we will derive below, and they have been esti- 
mated in observations over a range of redshifts. In com- 
parisons to published data, we obtain the relevant clus- 
ter mass assuming an NFW profile with a concentration 
parameter c = 5. We find that the data suggest a non- 
evolving (or- = 0) entropy floor of Kq — 125/i -1 / 3 keV 
cm 2 with a fixed gas fraction f g = 0.8, which are rough! 
con sistent with the observations oflPonman et al.l (|199' 
and lLlovd-Davies. Ponman. fc Cannon! (|2000f ). 

In Figure [2j in the lower panel, we show 
the local M^qq — Lh n i s caling relation found by 
iReiprich fc Bohringerl (|2002f) at < z < 0.1, along with 
the predictions from the unmodified model (dot-dashed 
curve) and from our fiducial model (solid curve) at 
z = 0.05. In the upper panel of the same figure, we 
show the observational data of Maugh an et alj (2005) 
for the same scaling relation over the redshift range 
0.7 < z < 0.9, compared to our predictions for both 
the unmodified and the fiducial models at z = 0.8. In 
the model predictions in both panels, we adopt the 
same cosmology as used in the observational papers 
(different from our fiducial cosmology, discussed below). 
The plot shows good agreement between the data 
and our fiducial model over the full redshift range 
(and the unmodified self-similar model is clearly ruled 
out). The models also predict the ICM gas temper- 
ature, and as another check, in Figure [3] we compare 
these predictions to dat a on the local T — M n nn an d 
T — Lboi scalings from IReiprich fc Bohringerl (2002). 
Again, we find rough agreement between our predic- 
tions and the observed data, although the models 
slightly overpredict the temperatures of the lowest-mass 
clusters. We find that our models are also in good 
agreement with other determinations of both t he local 
mass-temperature dWhite. Jones, fc Formanl I1997L 
Girardi et al.1 H998I: iFinoguenov. Reiprich fc Bohringerl 



20011 : ISanderson et al.l l2003f) and temperature- 
luminosity ([Markeyichl 119981: lArnaud fc Evrardl 119991 : 
iFairlev et all 120001 : iNovicki. Sornig. fc Henrvi 12001 
relations. 

The predictions of the fiducial model were then com- 
pared to measurements of the SZ decrement. In FigureSJ 
we show our model predictions at four different redshifts 
between 0.1 < z < 0.4, for both the M500 — S arc ,u/ fu and 
Lboi — S a rc,vl fu scalings, where S arc , v /fu is the frequency 
independent integrated SZ flux in a circle with radius 1 
arc-second centered on the cluster. Al so shown in the 
figure a re the observational data from IMcCarthv et al.l 
(2003b) and references therein, for the full redshift range 
z ~ 0.1 — 0.4. Again, we find good agreement between 
our fiducial model and the observed scalings involving 
the SZ decrement - although the data is still sparse, and 
shows a large scatter. 

For the fiducial cosmology, we adopt 
(n m h 2 ,£l DE , w , w a ,cr s „VL b h 2 ,n s ) = 
(0.14,0.73,-1.0,0,0.7,0.024,1.0). With the excep- 
tion of erg, these values are taken from the best-fit 
ACDM model foun d in the first-year WMAP data 
([Spergel et al.1 l2003h . It has been shown previously 
that the local cluster mass-temperature relation implies 
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a low er value of as than CMB observations (|Seliakl 
2002). We arrive at the same conclusion, when we 
adjust as such that the model predictions agree with 
the number counts of X-ray clusters. For the flux 
limit of F x (0.5 - 2.0 keV) > 3 x KT 1 * erg s" 1 cm" 2 
of th e deepest ex i sting X-ray surveys (fVikhlini n et al.l 
119981 iGioia et all 120011: iRosati et all [2002J), which is 
also approximat ely the limiting: flux proposed for a new 
all-sky survey (jHaiman et a l. 2005), we find that we 
need to set as = 0.7 in order to be consistent with 
the observed cumulative number counts (about 5.5 
clusters/degree 2 for the above threshold; see Fig. [5|). 
This low value of as is, in fact, consistent with th e 
recent year-three WMAP results (|Spergel et al.ll2005D . 
and also with previous es timates based on the local 
mass -temperature relation (Scliak 2 0021: iPierpaqli et al.l 
120031) . galaxy velocity fields (IWillick fc Strauss! 11998ft. 
and weak lensing measurements (iBrown et al.l 120021 : 
iHoekstra et~aH 120021 : Uarvis et al.l |2003j). Our iducial 
cosmology is also consisten t with other cluste r count 
measurements using optical (|Bahcall et al.ll2003f) and X- 
ray selected samples (iBorgani et al.ll2001l : IShuecker et al.l 
[20031: iGladders et al.N200ffl ~ 

Finally, despite its rough success, we call attention to 
potentially interesting tension between our fiducial model 
predictions and observations. In particular, for both the 
local and high redshift scaling relations shown in Fig- 
ure [31 the cluster temperatures are somewhat overpre- 
dicted for the lowest mass clusters. Likewise, Figure [5] 
shows that our model log N — log S relation somewhat 
overpredicts the number of bright clusters. That the 
current cluster data reveals such tension highlights the 
potential of larger surveys in calibrating models of the 
ICM simultaneously with cosmology. In the present pa- 
per, however, our goal is to explore a methodology, rather 
than a specific cluster model. We made use of the exist- 
ing data only to guide us in choosing a reasonable set 
of fiducial parameters, and we postpone a more rigorous 
study of the viability of this model to future work. 



4. THE PHENOMENOLOGICAL MODEL 

Our aim in this paper is to quantify the merits of us- 
ing a model for the cluster, as opposed to directly self- 
calibrating the m ass-observable relation, as proposed in 
previ ous studies (jMaiumdar fc Mohrl 1200 4: Wan g et all 
2004). The later approach will be hereafter referred to 
as a "Phenomenological Model." To enable a compari- 
son between the two approaches, we here consider joint 
cosmological and "cluster model" constraints, adopt- 
ing simple parametric power-law scalings for the mass- 
observable relation to in clude the latter. We follow 
iMajumdar fc Mohrl (|2004h and IWang eFall (|2004f L and 
parametrize the X-ray luminosity as 



L X , M = 4ttF x 4 = A X M^ 15 E 2 (z)(l + z) 



(21) 



where dh is the luminosity distance in Mpc, L x is the 
bolometric X-ray luminosity in units of erg s , F x is 
the observed bolometric X-ray flux in units of erg s _1 
cm~ 2 , M2oo,i5 is the mass of the cluster within an over- 
density of 200/3 c in units of 10 15 M , with a normaliza- 
tion \ogA x = —3.56 (see §5.1 for details). We apply a 
K-correction to obtain the flux in the 0.5-2.0 keV band, 



L x (0.5 - 2.0keV) = L x . bol k{z), with 
rMi+z) 



k(z) = 



(22) 



J °°A v (T)dv ' 

where A„(T) is the iRavmond fc Smith! <|1977f) cooling 
function for a gas of metalicity Z = O.SZq and tem- 
perature T. The ICM gas temperature in the direct self- 
calibration case is needed only for this K-correction, and 
was esti mated using the observed mass-temperatu re re- 
lation of iFinoguenov. Reiprich fc Bohringerl (2001), 



Mi 80 = (l±0.2)xl0 15 /i" 1 A/ G 



T 



3/2 



f3 T Av 3 keV. 



(23) 



where -Miso is the cluster mass within an overdensity of 
180 relative to the background, (3t = 1-75 ± 0.25, and 
A„ is the virial overdensity. 
The SZ flux is parametrized similarly as a power-law, 



S v d\ = A sz M^, 15 E 2/3 (z)(l + z)^ 



(24) 



where cLa is the angular diameter distance in Mpc, and 
S v has units of mJy at a frequency of v = 150 GHz, with 
a normalization \o%A sz — 8.29 (see §5.2 for details). 

5. SIMULATING GALAXY CLUSTER SURVEYS 

To simulate the abundance of clust ers of galaxies, we 
begin with the standard formalism (|Press fc Schechterl 
[1921, 

"'" / ( M )^ dl0g<r " 1 (25) 



din A/ 



M dhxM 



where M is measured out to an spherical overdensity of 
180 relative to the background, p is the present day 
background matter density, f (M) is taken f rom t he ex- 
tended Press-Schechter fit of I Jenkins et al.l ( |200ll ). and 
a is the variance of the linear density field at redshift 
z. The po wer spectrum was calcul ated using the fitting- 
formulae of lEisenstein fc Hul (|1999l ). modified to include 
the effects of a time- varying dark energy parameter with 
w(z) ^wn+w^ (|Linde r! !2003fl . 

The quantity we will use below to derive cosmological 
constraints is the number of clusters within a bin cor- 
responding t o a ra nge of the observable 4'. We follow 
iLima fc Hul (|2005f ) and include a constant log-normal 
scatter in the cluster mass at a fixed observable, with 
standard deviation eri ogM i^. In the absence of scatter, 
we use our mass-observable relation - either the phe- 
nomenological or physical model - to map the flux bin 
< * < to a mass bin M l obs < M obs < M*£\ to 

compute the expected counts. In the presence of scat- 
ter, for a redshift bin centered at zj of width Az, the 
expected counts within VE^ < \& < is modified to 

Ni ( Zj ) = AVAz— / ^ [erfcfo) - erfc(x l+1/ 



where 



log M - log M 



obs 



(26) 
(27) 



and dV/dzdVL is the comoving volume element, and erfc 
is the complementary error function. 
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5.1. Survey Parameters 

Our mock X-ray survey has a flux detection threshold 
of F x (0.5 - 2.0 keV) = 3 x 1CT 14 erg cm 2 s _1 , compa- 
rable to the limit reached in existing deep surveys, as 
well as the threshold in the al l-sky survey propos ed to 
the Dark Energy Task Force (|Haiman et al.l [20051 ). In 
our fiducial model, this yields a cluster surface density 
of ~ 5.5 deg -2 . Including scatter in the mass-observable 
relation increases the total number of clusters by ~ 10% 
due to preferential upscattering of a larger population 
of less massive clusters, particularly at high redshift as 
the mass function steepens. We furthermore impose a 
strict luminosity floor of Lfi — 3 x 10 42 erg s _1 in the 
observed band, r oughly corresponding to the luminosity 
of small groups (|Arnaud fc Evrardlll999l ). below which 
both observations and our physical model become unre- 
liable. The logarithmic slope of the fiducial mass vs. X— 
ray luminosity relation is set to f3 x — 1.807 t o match the 
obser ved scaling, with no evolution j x = (|Wang et al.1 
2004), and a normalization log A,. = —3.56. This then 
also matches the cluster surface density predicted in the 
fiducial physical model, and facilitates a fair comparison 
of the two approaches. 

The mock SZ survey is modeled afte r upcoming obser- 
vations by the So uth Pole Telescope (jRuhl et al.l 120041 : 
IWang et al.ll2004h . with a limiting SZ flux of S v = 3.0 
mJy at v = 150 GHz. This flux threshold produces a 
yield that again matches the cluster counts of 5.5 deg -2 
in the mock X-ray survey, allowing an easy compari- 
son. The fiducial power-law parameters are chosen - as 
before - to match the local scalings with no evolution 
(Psz = 1-68 with 7 SZ = 0), and the normalization is cho- 
sen to agree with the number of clusters predicted by the 
physical model (logA sz = 8.29). 

In both surveys, and for both the physical and phe- 
nomenological model, we consider three cases: (1) a sin- 
gle flux bin extending from the detection threshold up 
to arbitrarily high fluxes, (2) 20 flux bins with no scat- 
ter, and (3) 20 flux bins with a log-normal scatter that 
has a fiducial value of <7i og M|* = 0.1 for both the X-ray 
and SZ mass-observable relation (|Lima & Hu 2005), but 
is assumed to be a free parameter. For each of these, we 
consider two different descriptions of the entropy injec- 
tion history. First, we fit a power-law evolution of the 
entropy floor of the form Kq\z) = K (z = 0)(1 + z) aK , 
with Kq and «k as free parameters. Second, we consider 
arbitrary evolution of the entropy floor Kq, in which Kq 
is fit independently in each of 40 redshift bins of width 
Az = 0.05. In both of these cases, the fiducial entropy 
floor is set to Kq — 125/i. 1 / 3 keV cm 2 with no evolution. 
For reference, in Figure [6] we show the redshift evolution 
of dN/dzdft and the of the mass M t h corresponding to 
the detection threshold in both the X-ray and SZ sur- 



veys, each for the case of no scatter and 0~\ o 



gM|* 



0.1. 



Finally, for the purpose of comparison, the sky coverage 
in both surveys i s taken to be AO = 4000 deg 2 , the size 
planned for SPT (jRuhl et al.ll200l . For surveys covering 
a different solid angle, the constraints we obtain below 
scale as oc AQ' 1 / 2 . 

5.2. CMB Anisotropies 

Our simulated CMB survey is modeled after a near 
to medium-term space based all-sky CMB survey, sim- 



ilar to the proposed Planck mission with bands at 100, 
143, and 217 GHz. We assume fractional sky coverage of 
fsky ~ 0.8, and perfect subtraction of foregrounds. The 
Ct coefficient s are calculated up to £ max = 2 x 10 3 using 
KINK FAST (ICorasaniti eFipl200l . a version of CMB- 
FAST (jSeljak fc Zaldarringalll996h mod ified to include a 

Wa ng et alj (j2004D for more 



time varying w. 
details. 



See, e.g., 



6. THE FISHER MATRIX FORECASTS 
6.1. Background 

With a fiducial cosmology and a prescription for sim- 
ulating galaxy surveys, one can compute lower limits on 
statistical errors achievable on model and cosmological 
parameters. In our self-calibrating approach, the cluster 
observations are used to simultaneously constrain both 
the cosmology and model parameters. 

Assuming a reasonably well-behaved likelihood func- 
tion C, the Fisher matrix is a quick method to forecast 
j oint parameter uncertainties in a multi-parameter fit 
(jTegmark. Taylor, fc Heavensiri997h . The Fisher matrix 
is defined as 

where p a and pp are free parameters. The inverse of 
the Fisher matrix gives the best attainable covariancc 
matrix C. Therefore, the constraints on any parameter 
ps, marginalized over all other parameters, is 



Ap 2 = Css - [F- 1 ] 



(29) 



We further define a degeneracy coefficient £g, which 
quantifies the degree to which degeneracies among the 
parameters increase the marginalized parameter error 
over the single-parameter error, 



£,s — CssFss- 



(30) 



The degeneracy coefficient £g ranges from 1 for a purely 
non-degenerate matrix, to oo for fully degenerate param- 
eters. 

The combined Fisher matrix for a series of N separate 
experiments is the sum of the individual Fisher matrices, 
i.e. 



N 



F, 



(31) 



which motivates the definition of a synergy coefficient Q , 
to quantify the fractional improvement in the marginal- 
ized error on a parameter ps after combining experiments 
versus adding the individual experiment errors in quadra- 
ture, as 



\ " 



Cs = Css QQW) 



(32) 



The synergy coefficient ranges from 1 for no improvement 
from combining the experiments beyond adding the indi- 
vidual constraints in quadrature, to if the degeneracies 
are fully broken, in the limiting case that the combined 
Fisher matrix delivers a perfect knowledge of ps- 
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6.2. Galaxy Cluster Counts 
The Fisher matrix for a mock cluster survey, , can 
be calculated from the number counts of clusters as a 
function of the observable \&. We define the counts Fisher 
matrix as 

v dN z dNi 1 



F 



af3 



dp a dpp Ni 



(33) 



where i is the bin index, and the error is purely Poisson 
counting erro r (ignoring the modest increase due to sam- 
ple variance, iHu fc Kravtsovi (|2003l ). This represents a 
sum of several Fisher matrices, one for each independent 
flux and/or redshift bin. 

Therefore, given redshift information, it is possible to 
derive cosmological constraints from a combination of the 
flux distribution and redshift evolution of cluster counts. 
We bin the data by flux - as before - and redshift, giving 
a Fisher matrix 



F 



a/3 



EE 



dNi(zj) dNi{z 3 



1 



dp a dpp Ni(zj) 



(34) 



where the sum is taken out to a limiting redshift of zt. 
The first nu n — 1 flux bins are spaced evenly in log-space 
from the limiting flux over five orders of magnitude, with 
a final bin extending out to oo. We vary all the cosmolog- 
ical and model parameters, with the exception of fl^h 2 
and n s , which are later included via CMB observations. 
Finally, surveys in different observables are not assumed 
to overlap. 

6.3. CMB Observations 

In addition to constraints from cluster counts, we cal- 
culated the Fisher matrix for observations of CMB tem- 
perature and polarization anisotropies. We assume a 
Planck-like satellite mission 2 including both the tem- 
perature T and E-mode polarization auto-correlation, 
as well as the TE cross-correlation, while neglecting B— 
mode polarizatio n. The full CMB Fisher matrix is the n 
calculated as in [Zaldarriag a. Spergel. &; Seliakl (|1997f) . 
Each frequency channel was assu med to provid e indep en- 
dent cosmological constraints. See lWang et al.l (|2004! ) for 
full details. 

7. RESULTS AND DISCUSSION 

In order to present our results more clearly, we will pro- 
ceed in pedagogical order, considering successively more 
realistic assumptions. 

7.1. Cosmology-Only Constraints 

In this subsection, we study the limiting case in which 
the cluster model parameters are known perfectly. We 
begin with the most naive experiment, using a single flux 
bin (i.e. only total cluster counts above the detection 
threshold, as a function of redshift), and ignoring any 
scatter in the mass-observable relations. The constraints 
on the cosmology along with their associated degeneracy 
coefficients for X-ray and SZ surveys, and the synergy 
coefficient for the combination are presented in Table [TJ 
The physical model breaks many cosmological degenera- 
cies, and therefore gives much better constraints, for the 
X-ray and SZ surveys individually as compared to the 

2 see www.rssd.esa.int/index.php?project=Planck. 



phenomenological model. This is particularly true of the 
dark energy equation of state parameters wo and w a , 
which improve by a factor of ~ 2 in the X-ray and ~ 4 in 
the SZ survey. The single parameter sensitivities - the di- 
agonal elements of the Fisher matrix - for wo and w a are 
comparable in both approaches (physical or phenomeno- 
logical model) . We find that the degeneracy breaking in 
the physical model is due to better absolute sensitivity 
to Q m h 2 . This is because changing £l m h 2 while hold- 
ing flbh 2 fixed (and assuming a flat universe) changes 
the mean universal baryon fraction. Accordingly, the 
value of the observables for a cluster of fixed mass will 
vary with £l m h 2 . As a result, changes in fl m h 2 manifest 
themselves strongly in the redshift distribution of clus- 
ter counts in a flux-limited survey, even at fixed flbh 2 , 
when the physical, rather than phenomenological model, 
is used to calculate the threshold mass. This difference 
should be robust, in the sense that any physical model 
of cluster structure is likely to predict fi TO /^-dependence 
of the X-ray flux and SZ decrement (over and above the 
Q m h 2 dependence contained in the Hubble parameter 
E(z) in the power-law relations in equations [211 and [24]) . 

While using the physical model significantly improves 
the degeneracy coefficients for cosmological parameters 
in both the X-ray and SZ survey, it also removes the 
synergy between the two, and degrades the synergy coef- 
ficient for their combination. For example, £ Wo ~ 0.60 for 
the phenomenological model, as compared to C, WQ ~ 1 for 
the physical model. We illustrate this effect by showing 
the marginalized constraints in the fl m h 2 — Q,de plane in 
Figure [7] The error ellipses shown for phenomenological 
model (right panel) are larger than for the physical model 
(left panel), and they are also more complementary, with 
major axes pointing in different directions. The overall 
constraints on the cosmology are, however, are still bet- 
ter when the physical models are used. 

We next present constraints using 20 flux bins for both 
the X-ray and SZ surveys in Table Using multiple 
bins allows additional information to be extracted from 
the shape of the cluster mass function. As a result, the 
degeneracy coefficients for all parameters are substan- 
tially reduced in both surveys. In fact, all parameters 
become nearly non-degenerate, and the synergy coeffi- 
cients in both surveys using either the phenomenologi- 
cal or physical model are close to ~ 1. Therefore, these 
constraints reflect the individual parameter sensitivities 
of the surveys. These sensitivities are roughly equal in 
the two approaches for all parameters, with the excep- 
tion of £l m h 2 , to which the physical model is a factor 
of ~ 6 more sensitive. This is not surprising, since the 
cosmological dependence of the physical model reduces 
to that of the phenomenological model in the limit of 
K ~ 0, expect for the fl m h 2 dependence. We therefore 
confirm our interpretation of constraints using 1 flux bin, 
and illustrates quantitatively the power of the model to 
resolve Q, m h 2 . This sensitivity, under more realistic as- 
sumptions, will help the physical model break more de- 
generacies than the phenomenological model. 

Constraints for more realistic surveys, in which scatter 
in the mass-observable relation will degrade the informa- 
tion that can be extracted from the mass-function shape, 
are summarized in Table[3] These use 20 flux bins in both 
the X-ray and SZ, but also include a log-normal scatter 
as a free parameter. The scatter smears out the shape 
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of the mass function, and increases the degeneracy coef- 
ficients for all parameters and in both surveys, to values 
comparable to those using a single flux bin and no scat- 
ter. Overall, the constraints are within a factor of two 
of the 1-bin, no-scatter case, indicating that scatter de- 
stroyed most of the information from the mass function 
shape. The phenomenological case still has higher over- 
all degeneracy coefficients, which makes the X-ray and 
SZ constraints more synergistic in combination. On the 
other hand, the absolute constraints for all parameters 
are better when the physical model is used. 

We also find that imposing priors on the scatter does 
not significantly improve constraint s on the cosmology . 
This contrasts with the results of iLima fc Hul {2005). 
They, however, use a fixed mass threshold. For a flux lim- 
ited survey, constraints on the cosmology degrade when 
a log-normal scatter is introduced as a free parameter, 
not primarily because of degeneracies among the scatter 
and the cosmology, but rather because the scatter it- 
self washes out the cosmological sensitivity of the mass- 
observable relation. To test this hypothesis, we consider 
constraints on the cosmology and scatter using a fixed 
mass threshold versus a fixed flux threshold, both with 
1 bin. We find that when a fixed mass threshold is used, 
fixing the scatter improves cosmological constraints by a 
factor of 2 in most cases, and up a factor of 4 in some. 
However, when a fixed flux threshold is used, the im- 
provement is typically no more than 30%. This implies 
that the mass-observable relation, which is flattened by 
the inclusion of scatter, is driving the constraints on the 
cosmology. 

In summary, for the idealized case when the cluster 
structure or mass-observable relation is assumed to have 
no uncertainty, we find that the physical model leads 
to better sensitivity to dark energy parameters in both 
the SZ and X-ray surveys individually, and also in the 
case when the two surveys are combined. We find that 
this improvement arises because the physical model pre- 
dicts an increased sensitivity of the observables to fl m h 2 , 
which removes degeneracies between Q m h 2 and dark en- 
ergy parameters. Interestingly, the degeneracies are im- 
proved so substantially in the individual surveys, that the 
combination of the two surveys is less synergistic than in 
the phenomenological case. 

7.2. Self- Calibrated Constraints 

In this subsection, we present self-calibrated con- 
straints, in which the parameters of the phenomenolog- 
ical or physical model, are assumed to be a-priori un- 
known. We will focus on the extent to which this de- 
grades cosmological constraints. It is worth noting, how- 
ever, that such a self-calibration approach can yield tight 
constraint not just on the cosmology, but also on the 
properties of the clusters themselves. 

In Table IU we show self-calibrated constraints in the 
simplest case of 1 flux bin with no scatter. The re- 
sults show that the phenomenological model parameters 
are extremely degenerate with the cosmology, suggest- 
ing their ability to introduce changes in dN/dz in ways 
that mimic the effects of changing the cosmology. This 
has been noted previously: for example, uncertainties 
in the normalization of the observed mass-temperature 
relation result in large unc e rtainties in the value of erg 
measurements (|Seliakl 120021: iPierpaoli et al.ll2003f) when 



this value is inferred from the cluster temperature func- 
tion. The severe degeneracies between cosmological and 
mass-observable relation parameters are manifest in the 
high values of the degeneracy parameters, especially for 
erg, with of £ CTs ~ 100 and 600 for the X-ray and SZ sur- 
veys respectively. When the X-ray and SZ surveys are 
combined, it is possible to break some of these degenera- 
cies, but for us, and for the dark energy parameters in 
particular, the synergy coefficients Q Wo ,w a ~ lj indicating 
that uncertainty in the slope, normalization, or evolution 
of the mass-observable relations in the X-ray or SZ can 
mimic the effects of variations in these cosmological pa- 
rameters. 

The above contrasts sharply with the results we obtain 
using a physical model - as the bottom half of Table [4] 
shows, there is only minor degradation of the cosmologi- 
cal parameter constraints due to the uncertainties in f g , 
Kq, and ax ■ In this case, the X-ray parameter degen- 
eracies on average increase only by a factor of ~ 2 — 3 as 
compared to the case of perfect knowledge of the model 
parameters, while the SZ degeneracies increase by no 
more than 40%. This demonstrates that the physical 
model parameters f g , Kq, and ax do not easily mimic 
the effects of a change in cosmology. The comparison 
with the phenomenological model is especially dramatic 
for the dark energy equation of state parameters in the 
SZ survey, whose constraints are improved by a factor of 
~ 10. In addition, because these parameters simultane- 
ously determine both the X-ray and SZ mass-observable 
relations, and because they enter into the expressions 
for the flux differently (through p 2 g for the X-ray and p g 
for the SZ flux), the two surveys are more complemen- 
tary than in the case when these parameters are fixed. 
Therefore, the synergy coefficients are smaller than their 
analogues in Table [TJ 

In Table O we extract information from the shape of 
the mass function, by using 20 flux bins, while still as- 
suming no scatter in the mass-observable relation. As 
before, we find that the shape of the mass function 
helps to break many of the parameter degeneracies, with 
£ ~ 1 — 2 for all parameters, and thus these constraints 
again are indicative of the single parameter sensitivities. 
In the middle and the bottom set of rows in this table, 
we consider two different parametrizations of the entropy 
injection history in the physical model: power-law, and 
arbitrary evolution of Kq with redshift. In both cases, 
the degeneracy coefficients for all parameters arc signif- 
icantly lower than those in the direct self-calibration 
models. They are somewhat higher for arbitrary Kq- 
evolution by a factor of ~ 1.5, indicating that relax- 
ing the assumption of a power-law evolution can better 
mimic the effects of dark energy. However, in this case 
of arbitrary Kq evolution, the X-ray and SZ constraints 
are more synergistic in combination. Therefore, when 
both the X-ray and SZ results are considered, the ab- 
solute parameter constraints are degraded by no more 
than ~ 20 — 30% even when Kq can evolve arbitrarily 
with redshift. We also see, as in the previous section, 
far greater sensitivity to f2 m /i 2 when using the physical 
rather than the phenomenological model. 

Finally, the constraints in the most realistic among our 
mock surveys, which includes 20 flux bins and a log- 
normal scatter as a free parameter, is presented in Ta- 
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ble[6j As we found in the the cosmology-only constraints 
above, the uncertainty in the scatter increases the de- 
generacies for all parameters, in both surveys, and for 
both the phenomenological and physical model results. 
For the phenomenological model, we recover the comple- 
mentarity of the two surveys. However, both cases of 
the physical model give smaller degeneracy coefficients 
and better absolute constraints for all parameters. In 
particular, using the physical model improves Wq and w a 
constraints by a factor of 2 relative to the phenomeno- 
logical model. We also note that allowing arbitrary Kq- 
evolution still does not significantly degrade constraints, 
even in the presence of uncertain scatter, relative to a 
power-law i^o _ Gvolution. Interestingly, we find that this 
is because for a physical model with arbitrary Kq evo- 
lution, the X-ray and SZ surveys are exceptionally com- 
plementary. In particular, the synergy coefficient for the 
dark energy parameters ( Wo ,w a < 0.7, indicating that 
constraints improve by more than 30% (over adding the 
constraints in quadrature) when the X-ray and SZ sur- 
veys are combined. This is illustrated visually in Fig- 
ure [SI where we show marginalized constraints in the 
wo — w a plane. The figure shows a significant reduction 
in the size of the error ellipse when the two surveys are 
combined. In particular, this reduction is much more 
significant in the physical model (left panel) than in the 
case of the phenomenological model (right panel). 

7.3. Complementary Observations 

Finally, we contrast the constraints we find from 
the mock cluster surveys to measurement of CMB 
anisotropies (by Planck). The constraints and degener- 
acy coefficients we find for Planck individually are listed 
in Table [7l The table shows that the forecasts for the 
combined X-ray and SZ cluster surveys, even when we 
allow for unknown scatter, and for arbitrary entropy- 
injection history, compares favorably with Planck, in- 
cluding fl m h 2 which is known to be measured to exquisite 
precision with the CMB. 

We also investigate the utility of CMB measurements 
in breaking the parameter degeneracies inherent in the 
X-ray and SZ cluster surveys. Constraints for the com- 
bination of Planck with the X-ray and SZ cluster sur- 
veys (including 20 bins and uncertain scatter, as in 
Table [6j are listed in Table [H The inclusion of the 
Planck data adds substantial sensitivity to the individ- 
ual parameters, which results in a significant increase in 
the degeneracy coefficients. We find that when Planck 
data is assumed to be available, the combination of X- 
ray and SZ data, using a phenomenological model, is 
not particularly synergistic. All the synergy coefficients 
((CMB + XR + SZ) > 0.75, with constraints on the 
dark energy equation of state unimproved relative to 
adding the three individual experiment errors in quadra- 
ture. This implies that CMB observations resolve the 
degeneracies that were previously broken by the combi- 
nation of X-ray and SZ data. This is also the case for 
a physical model with power-law Kq evolution. How- 
ever, for arbitrary Kq evolution, we still find significant 
improvement in constraints for the combination of X- 
ray and SZ data, even in the presence of CMB data, 
with ((CMB + XR + SZ) < 0.8 for all parameters, and 
((CMB + XR + SZ) ~ 0.7 for dark energy equation 
of state parameters. This degeneracy breaking is only 



slightly weaker than in the absence of CMB data. In- 
deed, though the addition of Planck data to the indi- 
vidual surveys yields noticeably improved constraints on 
dark energy, when is included in the combined surveys 
it helps little with the wq and uu a constraints, since their 
marginalized uncertainties are dominated by the degen- 
eracies that remain with the cluster model parameters 
(rather than among cosmological parameters). 

8. SENSITIVITY TO THE ADOPTED PHYSICAL MODEL 

While the above results are promising, one may worry 
that they rely on the specific choice for the underlying 
physical model. In particular, if real data is fit with 
the wrong model, the derived cosmological parameters 
may suffer from a bias larger than their uncertainty we 
inferred above. Such potential biases could be directly 
quantified by creating mock data using various "true" 
models, and obtaining the best-fit to these data using 
different, parameterized, "wrong" models. The magni- 
tude of the bias, in general, will depend strongly on the 
pair of models that are compared - and hence, ultimately, 
on which way the real data falls. A full exploration of 
these biases also requires going beyond the Fisher ma- 
trix approach adopted here, and will be postponed to a 
future paper. 

Nevertheless, one can use the Fisher approach to quan- 
tify the impact of wrong model assumptions, by adding 
new degrees of freedom, representing uncertainties about 
the underlying physical model, and evaluating the degra- 
dation in the statistical errors. This degradation can be 
regarded as the systematic error arising from the choice 
of the wrong model, among a restriced set of possible 
models. 

Specifically, here we allow the gas fraction to evolve 
as f g ~ (1 + z) Q/ «, and the entropy injection to vary 
with mass as Kq ~ M KM and radius as Kq ~ (1 + x) Kx . 
To maintain our rough agreement with observations, the 
fiducial model is chosen to match that used for the previ- 
ous analysis, with ctf g — km = k x = 0, and the three ad- 
ditional parameters included in our self-calibrated Fisher 
forecasts for a power-law evolving entropy injection with 
20 flux bins and a log-normal scatter. 

We then compute the constraints for cosmological pa- 
rameters as before, to see the extended to which they 
are degraded by the addition of these new degrees of 
freedom. We find that in the X-ray survey, the cosmo- 
logical parameter constraints degrade by about a factor 
of ~ 2, while in the SZ survey, the constraints degrade 
by no more than 15%. This behavior is to be expected, 
as the X-ray luminosity scales as ~ p g , as compared to 
the SZ luminosity that scales as p g , and is therefore more 
sensitive to changes in the physical model. Likewise, the 
model parameter constraints degrade by ~ 60% in both 
the X-ray and SZ surveys. When both surveys are com- 
bined, the cosmological parameter constraints degrade 
by no more than ~ 20% - < 10% for the dark energy 
equation of state - while the model parameter constraints 
degrade by ~ 50%. At the same time, we are able to con- 
strain ctf g ~ the evolution of the gas fraction - to within 
0.14 and 0.11 from the X-ray and SZ surveys individ- 
ually, and 0.06 from their combination; roughly similar 
to constraints on the evolution of the entropy injection. 
Likewise, we are able to constrain the values of the new 
parameters to within ~ 0.02 from the X-ray and SZ sur- 
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veys individually, and 0.01 from the combination. 

Therefore, we conclude that the greatest impact of in- 
troducing additional flexibility to the model is to degrade 
theoretical constraints on the model parameters them- 
selves. The cosmological constraints, particularly for the 
SZ and combined surveys, are more robust to relaxing 
the model assumptions. Using a model with three ex- 
tra degrees of freedom, even in the case of the X-ray 
survey, which is particularly sensitive to the model be- 
havior, the theoretical constraints using a physical model 
are substantially better than those using a phcnomcno- 
logical fit. While these results are encouraging, a more 
complete treatment, allowing for a wider range of (non- 
power-law) models, would be necessary to estimate pos- 
sible systematic errors. 

9. FUTURE WORK 

Though these results above are promising, we here note 
some important reservations. Firstly, the goal in our 
study is to motivate a methodology, rather than a spe- 
cific cluster structure model. We believe that the im- 
provement in the constraints in individual experiments 
arises because adopting a physical model introduces new 
cosmology sensitivity that is not present when the mass- 
observable relation is parametrized directly. For the pre- 
heated model we consider here, the dominant effect is the 
new fi m /i 2 -sensitivity, which we expect to hold generi- 
cally, regardless of the details of a cluster model. Like- 
wise, the improvement in the synergy between the SZ 
and X-ray surveys arises because the physical model re- 
lates the X-ray flux and the SZ decrement (one cannot 
be changed without changing the other). This, again, 
should be a generic feature of any model that directly 
models the ICM gas distribution. A caveat to this last 
conclusion is that the model may have uncertainties, such 
as gas clumping or a strong radial dependence of the en- 
tropy injection Ko(r), which will changed the X-ray flux 
relative to the SZ decrement, which will tend to reduce 
the synergy between the two surveys. 

Nevertheless, before fits to real data can be per- 
formed, there is significant work to be done in mod- 
eling cluster structure and evolution more accurately. 
Modified entropy models, though very successful at 
predicting mass-observable rel ations and their evolu- 
tion, are incorrect in detail dPratt fc Arnaudl 120051 : 
iPratt. Arnaud. fc Pointecouteaul l2006f ). We also find 
some tension between our model predictions and ob- 
served local scaling relations and number counts (see 
Figures [31 and[5]). Furthermore, in addition to having 
models that fit cluster structure data accurately, one 
has to also know the correct cosmological dependence 
of the observables in these models. In the future, we en- 
vision that three-dimensional simulations, including the 
relevant non-gravitational and gas physics, at least in 
parametrized manner, and also including asymmetries in 
cluster structure, as well as deviations from hydrostatic 
equilibrium, will provide such models. 
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In addition to an improved model of cluster structure, 
more realistic constraints require a detailed treatment of 
the systematics of any real survey. Our work represents 
a "best case scenario", in which the X-ray and SZ sur- 
veys are idealized in many ways. With the introduction 
of scatter, and using bin sizes far wider than any reason- 
able observational error on flux measurements, the inclu- 
sion of systematic uncertainties in flux measurements will 
not significantly change our overall results. On the other 
hand, completeness and selection issues will degrade fi- 
nal errors, and may effect the complementarity we have 
found. As large-scale cluster surveys start taking data, a 
more thorough understanding of the combination of the 
systematic and statistical uncertainties will be required 
to interpret our results. 

10. CONCLUSIONS 

Using a Fisher Matrix approach, we have demonstrated 
that the use of a physically motivated model of the ICM 
in mock large-scale cluster surveys gives significantly 
better constraints on cosmological and model parame- 
ters, and better synergy between SZ and X-ray surveys, 
than one can obtain by directly parametrizing the mass- 
observable relation. In particular, when both the cosmol- 
ogy and model parameters are included in the fit, for the 
simplest case of pure cluster counts as a function of red- 
shift (see Table 2]) , the physical model yields constraints 
on the dark energy equation of state that are 2 times 
tighter and 3 times less degenerate in the X-ray, and 10 
times tighter and 15 times less degenerate in the SZ than 
those using a phcnomcnological model. If the shape of 
the mass function and scatter in the mass-observable re- 
lation are included and the entropy floor is taken to be 
an arbitrary function of redshift (see Tabled, the dark 
energy parameter constraints are 20% times tighter and 
2 times less degenerate in the X-ray, and 2 times tighter 
and 2 times less degenerate in the SZ than those using a 
phenomenological model. In addition, these constraints 
are up to a factor of two tighter than those from simply 
adding the individual experiment errors in quadrature, 
relative to a minor 20% improvement from combining 
constraints using a phenomenological model. These re- 
sult suggest that parametrized physical models of cluster 
structure will be useful when extracting constraints on 
both cosmology, and cluster structure itself, from future 
large-scale SZ and X-ray cluster surveys. 
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Fig. 1. — The pressure, density, and temperature profiles for a cluster with a total (gas + DM) mass of 10 14 h 1 Mq (left panels) and 
10 15 h~ 1 M e (right panels). The profiles are computed at rcdshift z = 0, assuming that the intra— cluster gas is in hydrostatic equilibrium 
in an NFW halo. The three curves in each panel correspond to different amounts of preheating: Kq = (no preheating; dash-dotted 
curves), Kq = 50 (dashed curves), and 100 keV cm 2 (solid curves; close to our fiducial choice of 125/1 1 / 3 keV cm 2 ). 
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Fig. 2. — The local (bottom panel) and high— redshift (top panel) mass-luminosity relation between mass and X— ray luminosity for 
the self-similar model without preheating (dot-dashed curves) and the fiducial model, with a non-evolving entropy floor Kq = 125/]. 1 / 3 
keV cm (solid curve s). The data points show measurements in the range < z < 0.1 (Rciprich & Bohringcr 2002) and 0.7 < z < 0.9 
(Maughan ct al. 2005). The model predictions are computed at corresponding mean redshifts of (z) = 0.05 and (z) = 0.8. 
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Fig. 3. — The local mass— temperature (bottom panel) and temperature-luminosity (top panel) relations for the self— similar model 
without preheating (dot-dashed curves) and the fiducial mo del with a non-evolving entr opy floor Ko = f25/i 1 / 3 keV cm 2 (solid curves). 
The data points are from local measurements at < z < O.f (Rciprich & Bohringcr 2002), and the model predictions are computed at the 
corresponding mean redshift of (z) 
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Fig. 4. — Model predictions at four different redshifts in the range < z < 0.4 for scaling relations involving the SZ flux, in the fiducial 
model with a non-evolving entropy flo or of Ko = f25/i 1 / 3 keV cm 2 . The data points with errors, shown in three different redshift bins, 
are taken from [McCarthy ct al. (2003a) and references therein. 
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Fig. 5. — The cumulative X-ray source counts pr edicted in our fiducial model, as a function of soft X— ray flux in the 0.5-2.0 keV band, 
compared to observations by Vikhlinin ct al. (1998) (dashed curve) bracketed by corresponding lcr error bars (solid curves). The fiducial 
model predictions are shown as filled squares with Poisson error bars. At the flux limit of F x > 3 X 10 — 14 erg s — 1 used in our analysis, we 
find a source density of 5.5 deg - 2 , in agreement with observations, while at higher fluxes, the model somewhat overpredicts the number of 
clusters. 
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Fig. 6. — The lower panel shows the evolution of the limiting mass M t h corresponding to the flux detection threshold, and the upper 
panel shows the number counts of clusters above this mass, in the mock X— ray (solid curves) and SZ (dashed curves) surveys. The pairs of 
curves in the top panel correspond to the case of no scatter in the mass-observable relation (lower curve) and a log-normal scatter in the 
mass, at fixed observable, of <7i og M|* =0.1 (upper curve). 
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Fig. 7. — Parameter constraints in the Q, m h? — Q,de plane for direct self-calibration of the mass-observable relation (right panel) and 
using the physical model (left panel). The constraints have been marginalized over all other cosmological parameters. The figure assumes 
a single flux bin (the shape of the mass function is not used) and no scatter in the mass-observable relation (see Table [TJ. In both panels, 
we show constraints from the mock X-ray (solid ellipses) and SZ (dashed ellipses) surveys individually, and from the combination of the 
two surveys (dot-dashed ellipses). Note the overall improvement in the constraints when the physical model is used, arising from breaking 
the w — O m h 2 degeneracies within each survey. 
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Fig. 8. — Parameter constraints in the wo — w a plane for direct self-calibration of the mass-observable relation (right panel) and 
using the physical model (left panel). The constraints have been marginalized over all other cosmological parameters. The figure assumes 
ri5 in = 20 independent flux bins, and a log-normal scatter in the mass-observable relation that is aprori unknown (with the fiducial value 
°f CT logM|* = 0.1; see Table(6]|. In both panels, we show constraints from the mock X-ray (solid ellipses) and SZ (dashed ellipses) surveys 
individually, and from the combination of the two surveys (dot-dashed ellipses). Note the improved degeneracy breaking in the physical 
model, when data from the two surveys are combined. 
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TABLE 1 

COSMOLOGICAL CLUSTER CONSTRAINTS FOR n hin = 1 AND NO SCATTER 







Ap(XR) 




Ap(SZ) 






Ap(XR + SZ) 


C,{XR + SZ) 


Phcnomenological Model 


AU m h 2 


5.01 x 10" 


-a 


2.4 


8.98 x 10" 


_3 


4.9 


2.80 x 10" 


-3 


0.64 




AQ DE 


7.44 x 10" 


-3 


10 


9.40 x 10" 


-3 


11 


2.43 x 10" 


-3 


0.42 




Awo 


0.10 




9.0 


0.17 




17 


5.14 x 10" 


-2 


0.59 




Aw a 


0.31 




5.5 


0.57 




11 


0.21 




0.78 




Acr 8 


1.31 x 10" 


-2 


15 


1.31 x 10" 


-2 


19 


3.69 x 10" 


-3 


0.40 


Physical Model 


AQ m h 2 


2.67 x 10" 


-3 


1.6 


2.10 x 10" 


-3 


3.0 


1.40 x 10" 


-3 


0.85 




An DE 


4.68 x 10" 


-3 


4.0 


4.25 x 10" 


-3 


5.5 


3.04 x 10" 


-3 


0.97 




Awo 


4.90 x 10" 


-2 


2.9 


3.35 x 10" 


-2 


3.0 


2.64 x 10" 


-2 


0.96 




Aw a 


0.19 




1.8 


9.89 x 10" 


-2 


1.8 


8.34 x 10" 


-2 


0.96 




Acrg 


4.33 x 10" 


-3 


4.5 


3.98 x 10" 


-3 


5.6 


2.75 x 10" 


-3 


0.94 



TABLE 2 



COSMOLOGICAL 


Cluster Constraints 


FOR nun = 20 


and no 


SCATTER 








Ap{XR) 


£(XR) 


Ap(SZ) 


OSZ) 


Ap(XR + SZ) 


C(XR + SZ) 


Phcnomenological Model 


AQ m h 2 


7.08 x 10" 4 


1.1 


1.22 x 10" 3 


1.3 


6.06 x 10~ 4 


0.99 




An DE 


8.07 x 10- 4 


1.5 


1.15 x 10" 3 


1.8 


6.32 x 10~ 4 


0.96 




Awo 


6.85 x 10" 3 


1.2 


7.27 x 10~ 3 


1.3 


4.89 x 10~ 3 


0.98 




Aw a 


6.45 x 10~ 3 


1.0 


9.05 x 10~ 3 


1.2 


5.07 x 10~ 3 


0.97 




Ao-s 


1.31 x 10" 3 


1.6 


1.28 x 10" 3 


1.9 


8.49 x 10~ 4 


0.93 


Physical Model 


AQmh 2 


1.44 x 10" 4 


1.0 


2.01 x 10~ 4 


1.1 


1.17 x 10~ 4 


0.99 




AQ DE 


7.10 x 10~ 4 


1.1 


6.15 x 10~ 4 


1.3 


4.53 x 10~ 4 


0.98 




Aw 


9.21 x 10~ 3 


1.1 


1.01 x 10~ 2 


1.2 


6.61 x 10~ 3 


0.97 




Aw a 


1.22 x 10~ 2 


1.1 


9.55 x 10" 3 


1.1 


7.20 x 10~ 3 


0.96 




A as 


1.06 x 10 -3 


1.2 


1.06 x 10 -3 


1.5 


7.01 x 10~ 4 


0.93 



TABLE 3 

COSMOLOGICAL CLUSTER CONSTRAINTS FOR n bin = 20 AND CTi ogM i^ = 0.1 LOG-NORMAL SCATTER 







Ap(XR) 




Ap(SZ) 




Ap(XR + SZ) 


((XR+SZ) 


Phenomenological 


AQ m h 2 


1.94 x 10" 3 


2.3 


5.08 x 10" 3 


3.7 


1.31 x 10" 3 


0.72 




Afl DE 


4.23 x 10" 3 


6.5 


7.21 x 10" 3 


9.1 


2.07 x 10" 3 


0.57 




Awo 


6.21 x 10- 2 


6.0 


9.15 x 10~ 2 


10 


3.94 x 10~ 2 


0.77 




Aw a 


0.23 


4.2 


0.32 


6.7 


0.16 


0.86 




Acr 8 


7.26 x 10" 3 


9.0 


8.57 x 10" 3 


13 


2.83 x 10" 3 


0.51 




Acrjffl. 


1.22 x 10~ 3 








1.20 x 10" 3 






Acrsz 






5.73 x 10" 3 




3.11 x 10~ 3 




Physical Model 


AQ m h 2 


1.55 x 10" 3 


1.6 


1.51 x 10" 3 


2.3 


9.99 x 10~ 4 


0.92 




An DE 


2.75 x 10~ 3 


2.7 


3.81 x 10~ 3 


5.1 


2.11 x 10~ 3 


0.95 




Awo 


3.62 x 10~ 2 


2.5 


3.01 x 10~ 2 


3.0 


2.02 x 10~ 2 


0.87 




Aw a 


0.12 


1.5 


9.38 x 10~ 2 


1.8 


7.16 x 10~ 2 


0.96 




Aa$ 


8.01 x 10" 3 


8.7 


4.19 x 10" 3 


6.1 


2.39 x 10~ 3 


0.64 




Ao-jfi? 


1.61 x 10~ 2 








9.90 x 10" 3 






Ao-sz 






2.45 x 10~ 2 




1.4 
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TABLE 4 

Self-Calibrated Cluster Constraints for n^j. 



1 WITH NO SCATTER 







Ap(XR) 


£(XR) 


Ap(SZ) 






Ap(XR + SZ) 


C(XR + SZ) 


Phenomcnological Model 




2.73 x 10~ 2 


13 


0.15 




82 


1.18 x 10 -2 


0.44 






5.97 X 10 


85 


0.17 




206 


2.19 X 10 


0.39 




Auio 


0.21 


18 


0.50 




50 


0.15 


0.80 






0.32 


5. / 


1.1 




21 


0.27 


0.88 




A(T8 


8.31 x 10~ 2 


98 


0.43 




627 


5.51 x 10~ 2 


0.67 




A log A x 


0.26 










0.23 






A(i x 


0.21 










0.15 






A7 X 


0.64 










0.41 






AlogA sz 






1.3 






0.14 






Ap 3z 






0.83 






0.10 






A 7sz 






2.1 






0.53 




Physical Model 


Af2 m ?i 2 


4.01 x 10" 3 


2.3 


2.18 x 10" 


3 


3.1 


1.66 x 10" 3 


0.87 


K (z) = K (l + z) a '< 


An DB 


1.37 x 10~ 2 


11 


4.65 x 10~ 


3 


6.1 


3.63 x 10~ 3 


0.82 




Au>o 


0.12 


7.0 


3.39 x 10~ 


2 


3.1 


2.84 x 10- 2 


0.87 




Aw a 


0.19 


1.8 


0.10 




1.8 


8.54 x 10- 2 


0.96 




Aa 8 


2.01 x 10~ 2 


20 


5.77 x 10~ 


3 


8.1 


4.28 x 10~ 3 


0.77 




A/ 9 


7.07 x 10- 2 




3.76 x 10- 


2 




2.25 x 10- 2 






Ai^o 


18 




6.2 






3.0 








0.36 




0.10 






6.03 x 10- 2 





TABLE 5 

Self-Calibrated Cluster Constraints for nu n = 20 with no scatter. 







Ap(XR) 


Z(XR) 


Ap(SZ) 


H{SZ) 


Ap(XR + SZ) 


C(XR + SZ) 


Phenomcnological Model 


AQmh 2 


7.36 x 10" 4 


1.2 


1.78 x 10- 3 


1.9 


6.53 x 10- 4 


0.96 




AQ DE 


8.64 x 10~ 4 


1.6 


1.22 x 10~ 3 


1.9 


6.74 x 10- 4 


0.96 




Awo 


7.19 x 10- 3 


1.3 


8.03 x 10~ 3 


1.4 


5.28 x 10- 3 


0.99 




Aw a 


6.53 x 10- 3 


1.0 


9.31 x 10~ 3 


1.2 


5.16 x 10- 3 


0.96 






1.76 x 10- 3 


2.2 


2.13 x 10- 3 


3.2 


1.30 x 10- 3 


0.96 




AlOgAo; 


1.89 x 10- 3 








1.74 x 10- 3 






Af3 x 


3.42 x 10- 3 








2.95 x 10- 3 






A 7x 


9.23 x 10- 4 








9.01 x 10- 4 






AlogA sz 






4.08 x 10- 3 




3.20 x 10- 3 






A/3 SZ 






5.69 x 10- 3 




4.39 x 10- 3 






A 7sz 






5.07 x 10" 4 




4.81 x 10- 4 




Physical Model 


AQmh 2 


1.48 x 10" 4 


1.1 


2.08 x 10- 4 


1.1 


1.17 x 10- 4 


0.98 


K Q (z) = K (1 + z)°>k 


AQ DE 


7.40 x 10- 4 


1.1 


6.16 x 10- 4 


1.3 


4.59 x 10- 4 


0.97 




Awo 


1.00 x 10- 2 


1.2 


1.05 x 10~ 2 


1.3 


6.87 x 10- 3 


0.95 




Aw a 


1.23 x 10- 2 


1.1 


1.05 x 10- 2 


1.2 


7.46 x 10- 3 


0.93 




Acr 8 


1.88 x 10- 3 


2.0 


1.21 x 10- 3 


1.7 


8.63 x 10- 4 


0.85 




A/ 9 


4.18 x 10- 3 




3.60 x 10- 3 




2.48 x 10- 3 






AK 


0.74 




0.91 




0.52 






Aa K 


2.19 x 10" 2 




1.64 x 10" 2 




1.15 x 10- 2 




Physical Model 


AQ m h 2 


1.95 x 10" 4 


1.4 


2.53 x 10" 4 


1.3 


1.32 x 10- 4 


0.86 


Arbitrary Kq(z) 


AQ DE 


8.34 x 10- 4 


1.3 


8.69 x 10~ 4 


1.9 


4.90 x 10~ 4 


0.81 




Awo 


1.15 x 10- 2 


1.4 


1.37 x 10~ 2 


1.7 


8.08 x 10- 3 


0.91 




Aw a 


1.62 x 10- 2 


1.4 


1.48 x 10~ 2 


1.6 


8.68 x 10- 3 


0.79 




Act 8 


2.08 x 10- 3 


2.3 


1.62 x 10~ 3 


2.3 


9.82 x 10- 4 


0.77 




A/ 9 


4.48 x 10- 3 




4.10 x 10~ 3 




2.70 x 10- 3 






< AK > 


1.3 




1.2 




0.68 





TABLE 6 

Self-Calibrated Cluster Constraints for n hin = 20 and ui og M|* = 0.1 log-Normal scatter. 







Ap(XR) 




Ap(SZ) 




Ap(XR + SZ) 


C(XR+ SZ) 


Phenomcnological 


AQ m h 2 


5.93 x 10~ 3 


6.9 


2.89 x 10~ 2 


21 


2.58 x 10" 3 


0.44 




AQ DE 


1.43 x 10~ 2 


22 


1.32 x 10~ 2 


16 


5.66 x 10~ 3 


0.58 




Awo 


7.88 x 10~ 2 


7.7 


0.11 


11 


5.55 x 10~ 2 


0.88 




Aw a 


0.27 


4.9 


0.34 


7.1 


0.19 


0.91 




Acrg 


1.40 x 10~ 2 


17 


1.56 x 10~ 2 


23 


8.85 x lO- 3 


0.85 




A log A x 


2.56 x 10~ 2 








1.72 x 10~ 2 






A/3 X 


2.09 x 10~ 2 








1.22 x 10~ 2 






A-y x 


9.67 x 10~ 2 








5.80 x 10~ 2 






Aa XR 


1.24 x 10~ 3 








1.23 x 10~ 3 






AlogA sz 






3.20 x 10~ 2 




2.12 x 10~ 2 






A/3 SZ 






4.77 x 10~ 2 




1.36 x 10~ 2 






A 7S 2 






3.32 x 10~ 2 




2.80 x 10~ 2 






Aa sz 






6.25 x 10" 3 




5.71 x 10" 3 




Physical Model 


AQ m h 2 


1.66 x 10" 3 


1.7 


1.55 x 10" 3 


2.4 


1.05 x 10" 3 


0.92 


K (z) = K (l + z) a K 


An DE 


2.86 x 10~ 3 


2.8 


4.11 x 10~ 3 


5.5 


2.22 x 10~ 3 


0.95 




Awq 


5.18 x 10~ 2 


3.5 


3.03 x 10~ 2 


3.0 


2.16 x 10~ 2 


0.83 




Aw a 


0.13 


1.5 


9.60 x 10~ 2 


1.9 


7.30 x 10~ 2 


0.95 




Ao-g 


9.40 x 10~ 3 


10 


5.31 x 10~ 3 


7.8 


3.33 x 10" 3 


0.72 




Aa X R 


6.47 x 10~ 3 








4.89 x 10~ 3 






Aasz 






7.55 x 10~ 3 




6.30 x 10~ 3 






A/ 9 


2.09 x 10~ 2 




2.67 x 10~ 2 




1.10 x 10~ 2 






AK 


5.0 




3.8 




2.2 






Aa K 


8. (9 X 10 




7 in — 2 
/.DO X 10 




4.3/ X 10 




Physical Model 


AQ m h 2 


2.12 X 10" 3 


2.1 


2.17 x 10" 3 


3.3 


1.19 X 10" 3 


0.79 


Arbitrary Kq(z) 


An DE 


3.66 X 10~ 3 


3.6 


7.11 x 10- 3 


9.6 


2.56 x 10~ 3 


0.79 




A«!o 


7.73 x 10~ 2 


5.2 


6.05 x 10~ 2 


5.9 


2.61 x 10~ 2 


0.55 




Aw a 


0.23 


2.7 


0.17 


3.3 


9.25 x 10~ 2 


0.68 




Acrg 


1.45 x 10~ 2 


15 


7.85 x 10~ 3 


11 


3.62 x 10~ 3 


0.52 




AaxR 


7.79 x 10~ 3 








5.61 x 10" 3 






A<r sz 






9.77 x 10~ 3 




6.74 x 10" 3 






A/ 9 


2.48 x 10~ 2 




3.94 x 10~ 2 




1.23 x 10~ 2 






< A Ko > 


2.2 




2.8 




1.5 





TABLE 7 

Cosmological Constraints for CMB Observations. 



Ap(PLANCK) ((PLANCK) 



AQ m h 2 1.21 x 10 -3 23" 

A(l DE 3.46 x 10~ 2 72 

Aw 0.32 166 

Aw a 1.04 152 

Act 8 4.11 x 10~ 2 88 

Afl b h 2 1.36 x 10~ 4 

An s 3.54 x 10~ 3 



23 



TABLE 8 

Self-Calibrated Cluster+CMB Constraints for n bin = 20 and <ri ogM i^ = 0.1 log-Normal scatter. 







Ap(XR) 


i{XR) 


Ap(SZ) 


Z(SZ) 


Ap(XR + SZ) 


aXR + SZ) 


Power-Law Model 


AQ m h 2 


5.93 x 10" 4 


11 


5.34 x 10" 4 


10 


3.35 x 10~ 4 


0.89 




AU DE 


3.17 x 10~ 3 


8.2 


5.53 x 10~ 3 


13 


2.07 x 10~ 3 


0.76 




Awo 


6.07 x 10~ 2 


32 


8.67 x 10~ 2 


46 


4.77 x 10~ 2 


0.97 




Aw a 


0.20 


30 


0.25 


37 


0.16 


1.0 




Ao-g 


4.87 x 10~ 3 


12 


1.12 x 10~ 2 


29 


3.78 x 10- 3 


0.85 




An b h 2 


1.10 x 10~ 4 




1.11 x 10~ 4 




1.07 x 10~ 4 






An s 


2.56 x 10~ 3 




2.49 x 10~ 3 




2.30 x 10~ 3 






A log A x 


2.21 x 10~ 2 








1.07 x 10~ 4 






A/3 X 


9.97 x 10~ 3 








2.30 x 10~ 3 






A~f x 


4.26 x 10~ 2 








1.67 x 10~ 2 






Acr XJ j 


1.22 x 10~ 3 








1.22 x 10~ 3 






A\ogA 3Z 






2.71 x 10~ 2 




1.82 x 10~ 2 






A/3 SZ 






1.32 x 10~ 2 




1.22 x 10~ 2 






A"f sz 






2.50 x 10~ 2 




2.28 x 10~ 2 












5.00 x 10" 3 




4.50 x 10~ 3 




Physical Model 


AQ m h 2 


5.72 x 10 -4 


10 


3.79 x 10" 4 


7.2 


2.97 x 10 -4 


0.97 


K (z) = K (l + z) a K 


AQ, DE 


2.39 x 10~ 3 


5.5 


3.61 x 10~ 3 


9.0 


1.85 x 10~ 3 


0.93 




Awq 


3.18 x 10~ 2 


16 


2.85 x 10~ 2 


15 


1.98 x 10~ 2 


0.93 




Aw a 


0.12 


17 


8.44 x 10~ 2 


12 


6.67 x 10~ 2 


0.97 




Aa 8 


5.91 x 10~ 3 


14 


5.05 x 10~ 3 


13 


3.03 x 10~ 3 


0.79 




An b h 2 


1.13 x 10~ 4 




1.05 x 10~ 4 




1.03 x 10~ 4 






An s 


2.54 x 10~ 3 




2.33 x 10~ 3 




2.26 x 10~ 3 






AtT XJ j 


5.59 x 10~ 3 








4.74 x 10~ 3 






Actsz 






6.39 x 10~ 3 




5.26 x 10~ 3 






A/ 9 


1.37 x 10~ 2 




2.59 x 10~ 2 




1.04 x 10~ 2 






AK 


3.4 




3.7 




2.1 






Aa K 


5.02 x 10 -2 




7.50 x 10" 2 




3.88 x 10~ 2 




Physical Model 


AQ m h 2 


6.70 x 10" 4 


12 


4.96 x 10" 4 


9.4 


3.04 x 10~ 4 


0.80 


K (z) = K (l + z)°"< 


ACl DE 


2.73 x 10~ 3 


6.3 


6.46 x 10~ 3 


16 


2.08 x 10~ 3 


0.83 




Awo 


4.45 X 10 


23 


5.53 x 10~ 2 


29 


2.38 x 10 


0.69 




Aw a 


0.18 


25 


0.15 


22 


7.98 x 10~ 2 


0.70 




Acr 8 


7.49 x 10~ 3 


18 


7.44 x 10~ 3 


19 


3.37 x 10~ 3 


0.64 




An b h 2 


1.16 x 10~ 4 




1.08 x 10~ 4 




1.04 x 10~ 4 






An s 


2.67 x 10~ 3 




2.45 x 10~ 3 




2.27 x 10~ 3 






AtT XJ j 


5.3 








6.8 






Acsz 






4.1 




2.8 






A/ 9 


7.03 x 10~ 3 




8.33 x 10~ 3 




5.51 x 10~ 3 






AK 


1.62 x 10~ 2 




3.52 x 10~ 2 




5.81 x 10~ 3 






Aa K 


13 




8.3 




1.12 x 10~ 2 





